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The Green function method developed in Ref. [S. B. Chernyshuk and B. I. Lev, Phys. Rev. E 81, 
041707 (2010)] is used to describe elastic interactions between axially symmetric colloidal particles 
in the nematic cell in the presence of the external electric or magnetic field. General formulas for 
dipole-dipole, dipole-quadrupole and quadrupole-quadrupole interactions in the homeotropic and 
planar nematic cells with parallel and perpendicular field orientations are obtained. A set of new 
results has been predicted: 1) Deconfinement effect for dipole particles in the homeotropic nematic 
cell with negative dielectric anisotropy Ae < and perpendicular to the cell electric field, when 
electric field is approaching it's Frederiks threshold value E ^ Ec. This means cancellation of the 
confinement effect found in Ref. [M.Vilfan et al. Phys.Rev.Lett. 101, 237801, (2008)] for dipole 
particles near the Frederiks transition while it remains for quadrupole particles. 2) New effect of 
attraction and stabilization of the particles along the electric field parallel to the cell planes in the 
homeotropic nematic cell with Ae < . The minimun distance between two particles depends on 
the strength of the field and can be ordinary for . 3) Attraction and repulsion zones for all elastic 
interactions are changed dramatically under the action of the external field. 



I. INTRODUCTION 

Colloidal particles in nematic liquid crystals (NLC) 
have attracted a great amount of research interest over 
the last decade. Anisotropic properties of the host fluid 
- liquid crystal - give rise to a new class of colloidal 
anisotropic interactions that never occur in isotropic 
hosts. The world of anisotropic liquid crystal colloids 
is much more varied than that of isotropic liquids. 

Experimental study of anisotropic colloidal interac- 
tions in the bulk NLC has been made in These 
interactions result in different structures of colloidal par- 
ticles, such as linear chains along the director field, for 
particles with dipole symmetry of the director, P, 
and inclined chains with respect to the director for 
quadrupole particles [1]-[H. Colloidal particles sus- 
pended at the nematic-air interface form 2D hexagonal 
structures That structures can be photochemi- 

cally switched with help of laser light Quasi two- 

dimensional nematic colloids in thin nematic cells form 
a rich variety of 2D crystals by using laser tweezers ma- 
nipulations. There are 2D hexagonal quadrupole crys- 
tals [1^, anti- ferroelectric dipole type 2D crystals 
[1^ \ l4j and mixed 2D crystals [iHj sandwiched between 
cell walls. 2D colloidal crystals assembled from chiral col- 
loidal dimers in twisted nematic cell was found in fl6| . 
Long ranged elastic interactions between colloids have 
been experimentally found to be exponentially screened 
(confinement effect) in the nematic cell across distances 
compared to the cell thickness L [T7| . 

Experimental results are reproduced by using the Lan- 
dau - de Gennes free-energy numerical minimization ap- 
proach well as molecular dynamics [l^j- 

Theoretical understanding of the matter in the bulk 
NLC is based on the multipole expansion of the direc- 
tor field and has deep electrostatic analogies [2l|-[3l|. In 
spite of some differences in these approaches, only one of 



them |2l| gives exact analytical quantitative result which 
has been proven experimentally in unlimited nematic LC. 
Authors of measured directly the dipole-dipole inter- 
action of iron spherical particles in a magnetic field and 
found it to be in accordance with [2l| , within a few per- 
cents accuracy for a dipole term. Authors of and Q 
have measured experimentally dipole-dipole interaction 
and found it to be in accordance with [2l| within about 
10% accuracy. This allows to justify main assumptions 
of |2l| for spherical particles in infinite nematic liquid 
crystal. 

Recently a proposal was made (32l . [33j to expand the 
description of colloidal particles on the confined NLC, 
which may be a generalization of the method [2l| for 
the bulk nematic LC. Using that approach elastic in- 
teractions between colloidal particles were found in the 
nematic cell and near one wall with either planar or 
homeotropic boundary conditions. The proposed theory 
[33I I fits very well with experimental data for the confine- 
ment effect of elastic interaction between spheres in the 
homeotropic cell taken from [13] in the range 1-lOOOfcT. 

Experimentally infiuence of the external fields on the 
nematic colloids were performed in papers [s^-fssj. In 
[3^ bidirectional motion of colloidal particles was ob- 
served in nematic cell under the action of ac electric 
field. Electrically driven multiaxis rotations of colloids 
was found in [33 • Transformation of the 2D colloidal 
hexagonal lattice at the nematic interface into the chains 
under the action of the magnetic field was found in [sgI ] . 
Electrically driven transformations of the Saturn ring di- 
rector configuration around the spherical particle into 
point hedgehog was experimentally observed in [38| and 
theoretically predicted in [l^. Qualitative investigation 
of the magnetic field influence on the elastic quadrupole 
interaction was made in [39j . 

Nevertheless there were no effective theoretical ap- 
proaches for quantitative treatment of the elastic col- 
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loidal interactions in the presence of the electric (mag- 
netic) fields. In this paper we propose the Green function 
method which enables to describe quantitatively elastic 
colloidal interactions between axially symmetric particles 
confined in nematic cell under the action of the exter- 
nal electric (magnetic) fields. The current method is the 
generalization of the method proposed in [s^, HI] for de- 
scription of the colloids in confined NLC. 

The outline of the paper is the following: In Sec. II 
we present the general Green function method for the 
description of colloidal particles in confined NLC in the 
presence of the external electric (magnetic) field. As an 
illustration of the method we find elastic dipole-dipole, 
dipole-quadrupole and quadrupole-quadrupole interac- 
tions in the unlimited bulk NLC in the presence of the 
electric field. We draw an analogy with previous results 
obtained in [s^ for planar nematic cell and propose a cor- 
respondence between electric field in unlimited NLC and 
nematic cell without any field. From this correspondence 
we propose a criterion for stability of the Saturn ring vs 
point hedgehog configurations around spherical particle 
in the nematic cell. Sec. Ill presents the energy of elastic 
interactions in the homeotropic cell placed in the external 
field. We consider cases of positive Ae > and negative 
Ae < NLC dielectric anisotropy and find that they are 
different in essence. The Ae < case lead to the decon- 
finement effect for dipole-dipole interaction when electric 
field approaches the Prederiks threshold value E ^ Ef. 
Dcconfincmcnt effect is the cancellation of the confine- 



ment effect which found in [17[ for quadrupolc interac- 
tion. In this section Sec. Ill we find as well new effect of 
attraction and stabilization of the particles along the elec- 
tric field E parallel to the cell planes in the homeotropic 
nematic cell for Ae < case. 

In Sec. IV we present the energy of elastic interactions 
in the planar nematic cell placed in the external field. 
Here three possible orientations of the field are considered 
(E||x, y, z axis). Significant changes of attraction and 
repulsion zones takes place in this case. And Sec. V 
represents conclusions. 



II. EFFECTIVE FUNCTIONAL FOR THE 
COLLOIDAL NEMATIC CELL PLACED IN AN 
EXTERNAL FIELD 

Here we will describe briefly the method developed in 
[3^ [33! and extend it to the case of presence of an exter- 
nal electric or magnetic field. Consider axially symmetric 
particle of the micron or sub- micron size which may carry 
topological defects such as hyperbolic hedgehog, discli- 
nation ring or boojums. In the absence of the particle 
ground, the non-deformed state of NLC is the orientation 
of the director n||z,n = (0,0,1). The immersed particle 
induces deformations of the director in the perpendicular 
directions = x,y. In close vicinity to the particle, 

strong deformations and topological defects arise, but be- 
yond them, deformations become small. In the paper (23 | 
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FIG. L (a)Homeotropic nematic cell of thickness L. (b) Pla- 
nar nematic cell of thickness L 



this area with strong deformation and defects was called 
the coat. Beyond the coat, the bulk energy of deforma- 
tion may be approximately written in the form: 



bulk 



K 

— I d^'xiyn^i 



with Euler-Lagrange equations of Laplace type: 



Ar 



(1) 



(2) 



Then director field outside the coat in the infinite LC 



has the form n^(r) = -I- 3c- 



3c 



with p and c being dipole and quadrupole elastic mo- 
ments (we use another notation for c with respect to the 
c in [2l|, so that our c = |c ). Actually p and c are 
unknown quantities. They can be found only as asymp- 
totics from exact solutions or from variational ansatzes. 
It was found in [2ll | that p — aa^ , c — —j3a^ with a being 
the particle radius and for instance a = 2.04, /3 = 0.72 
for hyperbolic hedgehog configuration [2l|. Experiment 
gives a = 2.05, /3 = 0.2 ± 0.1. Consider that we have 
found p and c and they are fixed. Then we can formulate 
an effective theory which describes particles and interac- 
tions between them, remarkably well. In order to find 
energy of the system: particle(s) -f- LC it is necessary to 
introduce some effective free energy functional F^jf so 
that it's Euler-La gran ge equations would have the above 
solutions. In the |2l| it was found that in the one con- 
stant approximation with Frank constant K, the effective 
functional has the form: 



F, 



eff 



K 



d^x 



(Vn^)2 



47rP(x)9^n^ ~ 47rC(x)9,a^i 
(3) 
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which brings Eulcr-Lagrangc equations: 

An^ = in [a^P(x) - 9,9^C(x)] (4) 

where P(x) and C(x) are dipolc and quadrupole mo- 
ment densities, fj. = x,y and repeated /i means sum- 
mation on X and y like = dxTi^ + dyUy. For 
the infinite space the solution has the known form: 
n^(x) = Jd^^'j^[-d'^P{^') + d'^d'M^')]. If we 



consider P(x) = pS{x) and C(x) = c6{x) this really 
brings n,(r) = + 3c^,ny{r) = + 3cf|. This 
means that efTective functional ([3]) correctly describes 
the interaction between particle and LC via linear terms 
-47rP(x)(9^n^ - 47rC(x)(929^7i^. 

When the external electric or magnetic field is applied 
we need to add one more term to the free energy ([3]) to 
take into account interaction between field and ncmatic 
LC. Thus, the effective free energy for the system under 
investigation is: 



I 

Ff}f j d\ I - y (en)2 - 47rP(x)a^n^ + 47r9, C(x) 9^ n^j (5) 



where k"^ = {AttK)~'^ AeE"^ if an electric field is applied 
and fc^ = AxH^ in the case of a magnetic field; 
Ae = e|| — ej_ and = X|| = X± Eire dielectric and 
magnetic anisotropies of the NLC respectively; e||E, H 
is the unit vector in the field direction. Hereafter we 
will consider only electric field for convenience. But all 
results can be applied as well to the magnetic field case 
by the substitution -^i?^ — )■ AxH^- In this paper we 
will consider only cases when the electric field is parallel 
to the one of the any coordinate axis {x,y,z). When the 
electric field is parallel to the z axis the Euler-Lagrange 
equations are: 

An^ = ^n^ + 47r [9^P(x) - 9,9^C(x)] (6) 

When the electric field is applied along one p = a; or 
p = y direction, then Euler-Lagrange equations are: 

An^ = -fc2<5p^7V + 47r [9^P(x) - 9,9^C(x)] (7) 

So in general case we have two different Euler-Lagrange 
equations for and Uy. 

In the case of confined nematic LC with the bound- 
ary conditions n^(s) = on the surfaces E (Dirichlet 
boundary conditions) the solution of EL equations has 
the form: 

n^(x)= / dVG,,(x,x')[-a;p(x') + a;a;c(x')] (s) 

Jv 

where is the Green function for ri^j . In the case of field 
direction e\\z we have only one Green function G = 
which satisfies equation (Ax — fc^)G(x, x') = — 47r(5(x— x') 
for X, x' G V (V is the volume of the bulk NLC) and 
G(x, s) = for any s of the bounding surfaces S. 

In the case of e\\p, p = x oi p = y field direction we 
have two different Green functions G^ which satisfy equa- 
tions (Ax-l-fc^(5p^)G^(x, x') = -47r(5(x-x') for x,x' g V 
(V is the volume of the bulk NLC) and G^(x, s) = for 
any s of the bounding surfaces S. 



The mathematical symmetry property Gp(x, x') = 
G^(x', x) can be proved for the Green functions satisfying 
the Dirichtle boundary conditions by means of Green's 
theorem [42| . 

Consider N particles in the confined NLC, so that 
-PW = Y.iPi^^^-'^i) and G(x) = Ci(5(x- x^). Then 
substitutfon ® into F^^^f"^ brings: F^^^f'^ = V^^^ + 

jj^nteracUon ^j^g^.^ jjself ^ ^ _ jjself ^ ^^^^ jjself 

interaction of the i-th particle with the bounding surfaces 
in the presence of the field U^'^^^ ~ U^'^l^ + ^dQ^ + ^qq ■ 
In the general case, the interaction of the particle with 
bounding surfaces, (self-energy part), takes the form: 

uif = -27r/^p2a^a;ij^(x„ x'j|x,=xi (9) 
t/j^'-'' = -2^A>c(9^a;9^i?^(x„x:;)+a;a^9,i/^(x„x^))|x.=x; 
u^Q^f = -27r/cc2a.9;a^a;if^(x„x:)|x.=x^ 

where ir^(x,x') = G^(x,x') — (we excluded the 

divergent part of self energy ). 

Interaction energy [/'"te™ctio« ^ ^^^^ ^mt_ jj^^.^ j/mt 
is the interaction energy between i and j particles: 
= Udd + UdQ + Uqq: 

Udd = -47rA>p'9^a;G^(x„ x^) (10) 



UdQ = -47rA-{pc'a^9;9^G^(x„x;) +p'c9;5,,a,G^(x„x;)} 



Uqq - -4^A'cc'9,a^9^a;G^(x„x;.) 

Here unprimed quantities are used for particle i and 
primed for particle j. Udd, UdQ, Uqq means dipole- 
dipole, dipole-quadrupole and quadrupole-quadrupole in- 
teractions, respectively. 
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Formulas © and (|10p represent general expressions for 
the self energy of one particle, (energy of interaction with 
the walls) , and interparticle elastic interactions in the ar- 
bitrary confined NLC with strong anchoring conditions 
np(s) = on the bounding surfaces and external elec- 
tric or magnetic field applied parallel to the one of the 
coordinate axes {x,y, z). 

It should be noted here that the difference between 
dielectric constants of the nematic host and suspended 
particles can give rise to an electrostatic interaction be- 
tween latter [s^ . For instance dipole-dipole induced 
electric interaction F^iectric can be roughly estimated us- 
ing the relation for two spheres immersed in an isotropic 
fluid H^: 

^electric = 1^7r£o£LC TZ J— (H) 

\e + 2eLcJ f 

where R is the radius of the droplets, Eq is the permittiv- 
ity of free space, and Elc is the dielectric constant of the 
particle, Elc is the averaged dielectric constant of the 



LC and r is the distance between spheres. In the whole 
paper we do not take into account this effect and focus 
only on the elastic interactions 'P elastic between colloidal 
particles arising from deformations of the director field. 

Let us consider for beginning briefly an infinite nematic 
liquid crystal with positive dielectric anisotropy Ae > 0. 
In this case the external field can be applied only along 
the no, i.e. along the z-axis, without great change of the 
initial director field no — (0, 0, 1) . The corresponding 
Euler-Lagrangc equations are 

An^ - k^n^ ^ 4^[a^P(x) - 9^9.C(x)] (12) 

where ^ = x,y ;k > and ?i^(x — > oo) = 0. Hence the 
Green functions are 

-fc|x-x'| 

G^(x,x') = ^ ^ (13) 

|x-x'| 

Then using PH)) it is easy to find that all interactions 
between particles are exponentially screened in this case. 



Udd pp' 



{(1-3 cos^ 9){1 + kr) + k^r^ sin^ 9} e"^'' 



(14) 



- {pc' - cp'f-^ { (15 cos2 - 9)(1 + fcr) + k^^G cos" - 4) - fcV^ sin^ 6] e"'^'- 



A-kK 



(15) 



f^QQ _ £^1(9 _ 90 cos^ 61 + 105 cos'' e){l + kr) + k'^r'^ (4 - 39 cos^ -)- 45 cos^ i 



+ fcV^ (1 - 9 cos^ 6* -t- 10 cos"* e) - k'^r'^ sin^ 6 cos^ eije"'^''' (16) 



where r = |x — x'| is the distance between particles and 
6 is the angle between r and z-axis. Such a screening 
of quadrupole-quadrupole interaction was qualitatively 
predicted in [39|. Maps of the attraction and repulsion 
zones of all interactions are presented on Fig [21 Fig (3] and 
Figgl 

These figures show big qualitative resemblance with 
maps of elastic interaction [Sj] (see as well thin black lines 
on FiglHl Figini FiglTT] for zero field) between colloidal 
particles in the center of the planar nematic cell with 
thickness L. So there is an analogy between infiuence of 
the electric (magnetic) field in unlimited NLC and con- 
finement of NLC with planar nematic cell. This resem- 
blance gives us a possibility to estimate correspondence 
between infiuence of the electric(magnetic) field in unlim- 
ited NLC and confined NLC without external field. For 
instance Fig|2] shows that deviation from the linear cone 
for dipole-dipole interaction and appearance of signifi- 
cant screening effects take place at distances k'^r'l « 2, 
Tc = ^/2/kE- For bigger distances screening effects ap- 
pear. As well Fig. 21 of [33| (and thin black lines on FigjS] 



a. as well) show that the same deviation takes place in 
the planar nematic cell with thickness L at the distance 
rl/L'^ « 1/2 so that we can formulate the correspon- 
dence between influence of the electric(magnetic) field in 
unlimited NLC and nematic cell without field in the form: 



(17) 



where S^e = ^/kE,^H = 1/fc/f arc the electric(magnetic) 
coherence lengths. This means that elastic proper- 
ties of the unlimited NLC in the presence of the elec- 
tric(magnetic) field are qualitatively the same as of the 
NLC confined with nematic cell with thickness L sa 
2£^E ~ Physically this means that half of the cell is 
approximately the same as the electric(magnetic) coher- 
ence length. 

H. Stark as well emphasized an an alog y between con- 
fined geometries and magnetic field in [4l[ and gave qual- 
itative estimates for transition from dipole to Saturn ring 
configuration in nematic cell with planar anchoring ob- 
served experimentally in [s^l . Using qualitative specula- 
tions he came to the estimate ~ L/2 — a where a is 
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FIG. 2. (Color online) Map of the attraction and repulsion 
zones for dipole-dipole interaction (|14|l between two parti- 
cles in the infinite NLC in the presence of the electric field 
E||z, Ae > as the function of the dimensionless distance kr, 




FIG. 3. (Color online) Map of the attraction and repulsion 
zones for dipole-quadrupole interaction (|15p between two par- 
ticles in the infinite NLC in the presence of the electric field 
E||2, Ae > as the function of the dimensionless distance kr, 
kE = E^Ae/4irK 



kz 




FIG. 4. (Color online) Map of the attraction and repulsion 
zones for dipole-dipole interaction (|16|l between two parti- 
cles in the infinite NLC in the presence of the electric field 
E\\z, Ae > as the function of the dimensionless distance kr, 
kE = Ey^Ae/'inK 

the radius of the particle that slightly differs from p7|) . 

In it was found that Saturn ring configuration in 
the presence of a magnetic field in unlimited NLC be- 
comes stable at a/£_H ~ 0.33. Then using analogy ([T7]) 
we can approximately estimate that Saturn ring configu- 
ration becomes stable in the nematic cell at the particle 
radius more than as'. 

as «0.16L (18) 

In the experiment [37| authors reported observation of 
the Saturn ring configuration around glass spheres with 
radii 20 and 50 fini confined in the center of planar ne- 
matic cell with thickness L = 120/^m. Then (fT8)) gives 
the lower limit of the radius as 05 — 19/im that allows 
experimental value of 20 . It would be really nice to 
test the ratio (|18p with more accurate experimental or 
computer simulation test. 

III. INTERACTIONS IN THE HOMEOTROPIC 
CELL PLACED IN THE EXTERNAL FIELD 

Let us choose coordinate system with z-axis parallel 
to the undistortcd director field ng and z = at the 
bottom wall. The x- and y-axis lie in the cell plane (see 
Figll^). There are two possible field orientations in this 
case. The first one is perpendicular to the cell walls, i.e. 
e = (0, 0, 1), and the second one is parallel to them, i.e. 
e = (1,0,0). 
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U/kT (a) U/kT (b) 




FIG. 5. (Color online) Log-log plots of the interaction energy in kT units as a function of the rescaled interparticle distance 
p/L. Here particle's radius a = 2.2 /im, cell thickness L = 7^m, K = 7pN, p = p' = 2.04a^, c — c' = Q.2a? . Dashed 
line 3 is the repulsion in the infinite NLC. Solid lines depict interactions in the homeotropic cell. Blue line 1 corresponds 
to the case Ae > and E = 0.99i?t. Green line 2 corresponds to the zero-field case E = G. Red line 4 corresponds 
to Ae < and E = 0.99i5t. (a) Dipole-dipole interaction. Red lines 4 -deconfinement effect; (b) Quadrupole-quadrupole 
interaction. Red lines 4 - no deconfinement effect. 



A. Field perpendicular to the walls 



And quadrupole-quadrupole interaction takes the form: 



Here we have E\\z, (en)^ = ~ 1 ^ 
Euler-Lagrange equations are: 



MM' 



The 



An^ - k^n^ = 4^[a^P(x) - d^d,C{^)] (19) 

and n^(z = 0) = n^(z = L) = 0, where ^ = x,y. Green 
functions can be easily derived from that one of the zero- 
field case [i^. They are: 



G^(x,x') = -£ 



L ^ ^ 

n—1 m—~oo 



sm — — X 



X sm ■ 



-Irn{XnP<)Krn{XnP>) (20) 



where Km are modified Bessel functions, tan(/3 ~ ^, 
tan ip' = Xn — \J + fc^. p< is the smaller of the 
p — \/ ^ and p' — \/ x''^ + y'^. Then (fTO|) brings 



dipole-dipole interaction in the homeotropic cell 

oo 



WttK ~ L 



riTTZ niTZ 
sm ^^sm^^ KaiXnP) (21) 



71=1 



L 



L 



here p ~ \/{y ~ v')"^ + (x — x')^ is the horizontal pro- 
jection of the distance between particles. It is easy to 
find that particles with z = z' repel when their dipole 
moments are parallel pp' > 0, and attract when pp' < 0. 
Similar dipole-quadrupole interaction is: 



UdQ_ 

16nK " L2 



1 °° 

j^'^ XlmrKoiXnp)^^ 



n=l 

rnrz' 



pc sm — — cos — — 



, mrz . niTZ 
cp cos sm 



(22) 



Uoo cc' v-^ , 9 9 9 nnz nnz' , 
I^^I^E^n^ '^^os — cos-r^o(^»^) (23) 

n—1 

When product cc' is positive quadrupole-quadrupole in- 
teraction between particles with the same ^-coordinates 
is repulsive, when cc' < interaction it is attractive. 

Application of ((2T|) for two spherical particles with 
z = z' = -J is shown in Figl5^. In the limit of small 
distances p <^ L interaction energy tends to the corre- 
sponding value in the unbounded NLC regardless of the 
field strength and nematic's anisotropy: Udd — ^ ^'"^pp _ 
If the cell is in the zero-field (see FiglS^, line 2) then for 
distances p > 1.2L exponential decay of interaction takes 
place [3^ . 

For further 



Let us take a look at A„ = 



L2 



P. 



analysis it is convenient to rewrite it as 




sgn(Ae) f E 
E't 



(24) 



where Et is the Fredericks threshold electric field, Et = 
T\l ^Ke\' '-''^^ ''^^ compare formulas ((2T|) - ((23)) with re- 
sults of [32I and find that an external field gives rise to a 
set of the effective cell thicknesses 



n^Ef 



n^Ef + sgn(Ae)£;2 



(25) 



and thereby affects the interaction between particles. For 
zero field £' = we have L'^ = L. It is clearly seen that 
when Ae > the external electric field decreases effective 



7 



thickness L^^ => (i? ^ oo) of the cell L and increases 
screening strength (see FiglS^, line 1). 

A completely different picture we can see in the ne- 
matic with negative anisotropy Ae < 0. In this case the 
maximum value of the applied electric field is the Freder- 
icks threshold value E = Et , after which Fredericks tran- 
sition takes place and initial ground state no = (0, 0, 1) 
becomes broken. So allowed values of the electric field are 
< E < Et (remember that we are bounded with the one 
constant approximation Ki ~ K2 = K3 = K and strong 
anchoring on the cell planes). So gradual increasing of the 
field strength below Et makes L'^ ^ 00 that is the De- 
confinement effect for dipole-dipole interaction between 
dipole particles in the center of the nematic cell near the 
Fredericks transition point (see FiglS^, line 4). Here at 
the distances 0.7L < p < 8L energy of dipolc-dipole in- 
teraction decreases even more slowly than in the infinite 
NLC. Actually this means that influence of the bound- 
ing surface planes disappears near the critical point Et 
for dipole-dipole interaction. Actually all L'^ for n > 2 
remain finite. It is interesting that this Deconfinement 
effect is absent for quadrupole-quadrupole interaction be- 
tween particles in the center of the cell (see Figl5]3, line 
4). Actually shows that the first nonzero term in 
this case is n = 2 , not n—1, and all for n > 2 
remain finite. So approaching the Fredericks transition 
only weakens the screening effect a little bit (FiglSjD). 



And quadrupole-quadrupole interaction is: 
^ 2^ n^TT^ cos — cos -J- (A„ + cos 29) 

n—1 

(29) 

In previous Sec lIII Al we have found that interaction be- 
tween particles with z = z' = ^ and positive products 
pp' > 0, cc' > is isotropic and repulsive throughout the 
cell plane. Expressions ^7} and show that external 
field, when it is applied parallel to the walls, breaks the 
symmetry of interaction and induces zones of attraction 
(see Fig|B^ and|B)D). Energy has a minimum at the points 
Pq of intersection between zones borders and coordinate 
axes (see Figint and [BJl) . Hence the distance po between 
the center and intersection points is the equilibrium dis- 
tance between particles. It depends on the field strength: 
the stronger the field, the smaller the po- Figure [7] shows 
that the particles suspended in a nematic with negative 
anisotropy can be can be arbitrary close. When Ae > 
the field strength must be less than Et to avoid the Fred- 
ericks transition, therefore, for dipoles we have po > 0.8 L 
and for quadrupoles po > 4.1 L. This effect of attraction 
and stabilization of the colloidal particles is new and have 
never been reported before. 

IV. INTERACTIONS IN THE PLANAR CELL 
PLACED IN AN EXTERNAL FIELD 



B. Field parallel to the walls 



Here we consider E||a;, (en)^ 



So that we have 



different Euler-Lagrangc equations for rix and n^: 
An, + k^n, = 47r[9,P(x) - d,d,C{x)] 

Any - 4^[ayP(x) - dyd,C{x)] 



(26) 



with boundary conditions n^{z = 0) = n^{z — L) — 0, 
p = x,y. To get the corresponding Green functions one 
can use ([^0]) with substitutions k'^ — — fc^ and fc^ = for 
Gx and Gy, respectively. Then ([TO)) brings dipole-dipole 
interaction: 

Udd pp' . nirz . mrz' . , , „ 
— — = — > sm — — sm — — [An + cos 29) (27) 
SttK L ^-^ L L 

n—1 

where 9 is the azimuthal angle between p and x. 

An = A^A'o(A„p) + PnK^iPnP) 
B„ = XlK2{XnP) - PnK2ipnP) 



with p = ^ and A„ 
interaction is: 



^^rf^ — fc2. Dipole-quadrupole 



^ 00 
T^-nY^riTT {A„ 



8ttK 



Bn COS 29) X 



X I pc sm — — cos — — 



, mrz . mrz 
p c cos sm 



(28) 



Coordinate system is depicted in FiglTJa. Here we have 
three possible field orientations. The first one is perpen- 
dicular to the ceU walls, i.e. e = (1,0,0). The second 
one is parallel to the walls and perpendicular to no, i.e. 
e = (0, 1,0). And finally, external field may be parallel 
to the ground director state no, i.e. e = (0,0, 1). 



A. Field perpendicular to the cell planes 

Here we consider 'E\\x, (en)^ ~ n^. The Eulcr- 
Lagrange equations are the same as (j26p . But have 
another boundary conditions (z is replaced with x), so 
that n^[x = Q) = n^(x = L) = p = x,y . Hence, the 
Green functions can be constructed from ((20|) replacing 
(x, y, z) — )• (y, z, cc) as well as fc^ — fc^ and fc^ = for x 
and y respectively, so that: 



.00 00 

4^ ^ e-(^-^')sin!^x 



G.(x,x')=^^ y: 



n—1 m— — 00 



L 



X sill ——Im{XnP<)Km{\nP>) 



00 00 



G.(x,x') = yE E 



n—1 m— — 00 



L 



e"""' ' sm — :— X 
m:x' 



{pnP<)Krn{p„P>) (30) 
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FIG. 6. (Color online) Map of the attraction and repulsion zones of interaction between two particles in the middle of the 
homeotropic cell under external field E = 0.99-Et. Field is parallel to the cell walls E||a::. Blue lines 1 correspond to the positive 
Ae, red lines 2 correspond to the negative Ae. Dashed lines depict zones in E — 2Et and Ae < 0. Sign "-" means attraction, 
"+" means repulsion, (a) Dipole-dipole interaction, pp' > 0. (b) Quadrupole-quadrupole interaction, cc' > 0. (c) Dimensionless 
energy of dipole-dipole interaction, (d) Dimensionless energy of quadrupole-quadrupole interaction. E = 2Et and Ae < in 
both graphs. 



where A„ 



of the p ~ y/ z'^ + y'^ and p' ~ \/ 



j^. p< is the smaller 



tant/j' = ^7. Then expressions ([TU)) give us all desired 
interactions in the planar cell. Dipole-dipole elastic in- 
teraction is given by: 



with F\ and F2 being: 



L /i„ . nixx , n-Kx' 



L 



L 



[Ko{p.nP) + K2{p-„p)] ■ 



n TT cos ■ 



mrx TiTTX 
■ cos — — - 



Ko{Xnp) 



F2 - ^ LV' sin — sin -— /C2(/i„p) (32) 



L L 



where p ~ — y'Y + (-^ — z'Y is the horizontal pro- 
IQnKpp' 2 jection of the distance between particles and 6 is the az- 

Udd = (^1 — F2 cos 6*) (31) imuthal angle between p and z. 
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FIG. 7. (Color online) Equilibrium distance between dipoles 
(blue lines 1) and quadrupoles (red lines 2) as a function of 
the field strength. 



Dipole-quadrupole interaction takes the form: 
16ttK 



where 



L4 



(pc' - cp') cos 6 (Ci + C2 cos^ e) (33) 



C2 



2F2 
p 



Quadrupole-quadrupole interaction is: 



where 



L5 



Di + D2 cos^ + D3 cos^ e) (34) 



P 



Do = L' - 



IOF2 



As = L' 



8F2 

/i2 



p 

5R 



2p 



PP 



Maps of the attraction and repulsion zones for interac- 



tion between particles with x 



are presented in 



FiglH In the limit of small distances p <C L all interac- 
tions are the same as in the unbounded nematic: Udd ^ 
^^(l-3cos2 0), C/dQ ^ 4H^(pc'-cy)(15cos20-9), 

f^QQ -> ^^^j^ (9 - 90 cos^ 6* 105 CDs'* 0) . When p > lOL 



zone's borders coincide with those in the cell under zero- 
field. At the intermediate distances scale field moves 
zone's borders closer to those in the unbounded NLC 
when Ae > (see Fig HI lines a). If the nematic has a 
negative Ae then field shrinks parabola-like zones and, in 
fact, enhances effects caused by the bounding walls (see 
FigEl lines b). 



B. Field parallel to the cell planes and 
perpendicular to the rubbing direction 

Let us consider the case with the field is parallel to the 
cell planes and perpendicular to the rubbing direction no, 
i.e. E||y, e = (0,1,0). 

In order to get corresponding Green functions we can 
simply make a replacement x O y in (j30p . One can eas- 
ily find that such replacement means a replacement of 

A„ = ~ with /i„ = ^ and vice versa in (j32[) . 

Formulas (|3T|) . ((33)) . (p4)) are valid in this case as well as 
in the previous one. And their application to the in- 



teraction between two particles in the center of the cell 
X = x' = ^ is presented in FiglHl It is seen that the 
external field does not affect the interaction at small dis- 
tances p ^ L. Though it changes maps of attraction 
and repulsion zones on far distances. Particularly, when 
the dielectric anisotropy is positive Ae > (see FiglHl 
lines a), the electric field moves zone's borders closer to 
their analogues - cones - in the infinite nematic LC. This 
means decreasing of the cell planes confinement effect as 
electric field tends to rotate the director from z direction 
to y direction. This approaching E ^ Et produces the 
Deconfinment effect in the planar cell similar to the de- 
confinement effect in the homeotropic cell described in 
Sec. Ill, a (see Fig. [T6| . And the Frederiks transition 
takes place at E = Et which breaks ground state no\\z. 

More interesting changes of the interaction maps take 
place for negative dielectric anisotropy Ae < 0. We have 
no limitations on the field strength E in this case. The 
electric field widens lateral tails of the dipole-dipole in- 
teraction with increasing of the strength value E (see 
Fig. (TU]). But there is the critical value Ecoi = 1.74£^t 

(Et = which a collapse of the lateral tales 

into dumbbell-shaped region occurs (compare Fig. [TT] 
and Fig. [T2). Particles repel inside the dumbbell-shaped 
regions and attract outside of them (see Fig. [TT|) . So 
that the attraction in the perpendicular to ng direction 
occurs a,t E > l.TAEt- Subsequent increasing of the 
field strength E compresses dumbbell-shaped regions and 
makes perpendicular attraction more significant. The 
same situation takes place for quadrupole-quadrupole in- 
teraction (see Figs. I13ll4ll5p . It is interesting that the 
same collapse and origin of the dumbbell-shaped regions 
occurs at the same field strength Ecoi = 1.74£'f as for 
dipole-dipole interaction. 
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(a) 
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FIG. 8. (Color online) Map of the attraction and repulsion zones of interaction between two particles in the middle of the planar 
cell with rubbing along z and thickness L placed in the perpendicular field E||a;, E — 0.99Et. Dashed lines are zone's borders 
in the unlimited NLC. Thick red lines a correspond to the positive Ae, thick blue lines b correspond to the negative Ae. Thin 
black lines are zone's borders in the zero-field. Sign "-" means attraction, "+" means repulsion, (a) Dipole-dipole interaction, 
pp' > 0. (b) Dipole-quadrupole interaction. The larger particle is in the center, pc' > cp' . (c) Quadrupole-quadrupole 
interaction, cc' > 0. 




FIG. 9. (Color online) Map of the attraction and repulsion zones of interaction between two particles in the middle of the planar 
cell with rubbing along z and thickness L placed in the electric field applied along the y-axis, E||y (see Figllja), E — 0.99Et. 
Dashed lines are zone's borders in the unlimited NLC. Thick red lines a correspond to the positive Ae, thick blue lines b 
correspond to the negative Ae. Thin black lines are zone's borders in the zero-field. Sign "-" means attraction, means 
repulsion, (a) Dipole-dipole interaction, pp' > 0. (b) Dipole-quadrupole interaction. The larger particle is in the center, 
pc' > cp' . (c) Quadrupole-quadrupole interaction, cc' > 0. 



C. Field parallel to the cell planes and parallel to 
the rubbing direction 



Finally, let's consider the case when the electric field E 
is parallel to the cell planes and parallel to the rubbing 
direction no, i.e. E||z, e = (0,0,1). In this case x and 
y directions are equivalent and we can find the Green 
function from (|30|1 as G = Gx ~ Gy with = = 

y^^^pH-"fc^. Then all interactions are again given by for- 
mulas dsn, (133), (El with fin = A„ = J^ + k^. Maps 



of the interaction zones are depicted on FigfTTl The elec- 
tric field shrinks zones for positive dielectric anisotropy 
Ae > and expands them to the unlimited nematic cones 
for negative Ae < 0. As in the Sec. Ill, A one can com- 
pare correspondent formulas with results of [32j and find 
that an external field gives rise to a set of the effective 
cell thicknesses 

^ ^^n^E^ +lgn{K^ ^^^^ 

and thereby affects the interaction between particles. For 
zero field £' = we have = L. It is clearly seen that 
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FIG. 10. (Color online) Dynamics of the attraction and repul- 
sion zones of dipole-dipole interaction between two particles 
in the middle of the planar cell with Ae < in the electric 
field applied along the i/-axis, E||y with increasing of the field 
strength E. The upper right quarter of the map is depicted. 
Dashed line is the zone's border in the unlimited NLC. Thin 
black line 1 corresponds to the E = Q, thick blue line 2 - 
E = Q.mEt, thick green line 3 - £ = l.&Et, thick black line 4 

- E = l.liEt, thick red line 5 - £ = l.TAEt, thick brown Une 6 

- E — 2Et. Attraction zone is above each line, repulsion zone 
is below each line, pp' > (see Fig llll Fig |12l as well). Col- 
lapse of the lateral tales into dumbbell-shaped region occurs 
at _E = 1.74Et. 



I rubbing 
— E 




FIG. 11. Map of the attraction and repulsion zones of dipole- 
dipole interaction between two particles in the middle of the 
planar cell with Ae < placed in the electric field applied 
along the j/-axis, E||i;, E = 1.73-Et. Dashed lines are zone's 
borders in the unlimited NLC. Sign "-" means attraction, "-|-" 
means repulsion, pp' > 0. 




FIG. 12. (Color online) Map of the attraction and repul- 
sion zones of dipole-dipole interaction between two parti- 
cles in the middle of the planar cell with Ae < placed in 
the electric field applied along the y-axis, E||y, E = l.lAEt- 
Dashed lines are zone's borders in the unlimited NLC. Col- 
lapse of the lateral tales into dumbbell-shaped region occurs 
at Ecoi ~ 1.74£'t (see Fig[T2]). Sign "-" means attraction, 
means repulsion, pp' > 0. 



when Ae > subsequent increasing of the field strength 
E will just make red zones on Fig ll7l ad lib narrow as 
it just decreases the effective cell thickness L„ and 
increases screening strength. 

Of course there is the limitation E < Et and Frederiks 
transition takes place when Ae < 0. In this case as well 
Lf^ — ^ oo as E ^ Et so tliat the deconfinement effect 
takes place before the Frederiks transition point. 

We see that maps of interactions on the Fig. [T7] are 
very resemble with Fig. ^ Fig. 21 So we can try to 
extend our speculations in Sec. II for the planar nematic 
cell with the presence of the field. As we have already 
mentioned in the Section II, the external field applied 
in the NLC with positive Ae > along the no induces 
the hedgehog transition to the Saturn-r ing , i.e. dipole 
type defect transforms into quadrupole [3 [13, HI] . At 
the same time the Ae < case leads to the deconfine- 
ment effect and the influence of the cell planes disappears 
near the critical point Et for dipole-dipole interaction. 
In some sense the confined NLC becomes similar to the 
unlimited nematic. But we know that in the unlimited 
nematic LC the hedgehog configuration is the most sta- 
ble [ll]. So there should be the point where Saturn ring 
configuration loses it's stability, i.e. has higher energy 
than hedgehog. In the formulas ([5T|) . ([55)1 . the main 
contribution has the term n = 1, so we may put forward 
a hypothesis that it's possible to replace L with L^^^ 
in the ([TS]) and to estimate the point of stability of the 
Saturn ring configuration in the planar nematic cell with 
external electric field: 
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FIG. 13. Map of the attraction and repulsion zones of 
quadrupole-quadrupole interaction between two particles in 
the middle of the planar cell with Ae < placed in the elec- 
tric field applied along the y-axis, E| E = 1.73Et. Dashed 
lines are zone's borders in the unlimited NLC. Sign "-" means 
attraction, "+" means repulsion, cc' > 0. 



FIG. 15. (Color online) Map of the attraction and repulsion 
zones of quadrupole-quadrupole interaction between two par- 
ticles in the middle of the planar cell with Ae < placed in the 
electric field applied along the y-axis, E| E — 2Et. Dashed 
lines are zone's borders in the unlimited NLC. Compression 
of the dumbbell-shaped regions occurs with increasing of the 
field strength for E > 1.74i5t (see Fig[T3}. Sign "-" means 
attraction, means repulsion, cc' > 0. 
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FIG. 14. (Color online) Map of the attraction and repulsion 
zones of quadrupole-quadrupole interaction between two par- 
ticles in the middle of the planar cell with Ae < placed in 
the electric field applied along the y-axis, E| E = 1.7 AEt- 
Dashed lines are zone's borders in the unlimited NLC. Col- 
lapse of the lateral tales into dumbbell-shaped region occurs 
aX E = 1.74Et (see Fig[T3j. Sign "-" means attraction, "+" 
means repulsion, cc' > 0. 



as « 0.16Ll^^ = 



0.16L 



l + sgn(Ae)( 



(36) 



For particles with radius a > as the Saturn ring has 
the less energy than hedgehog and vice versa. When 
Ae > the external field decreases effective cell width 
and strengthens Saturn ring stability. If the particle 
has radius a < 0.16L and has the hedgehog director 
configuration, the external field with value E > Ec = 

Et^Jl+J^^^, {Et = should make Saturn 

ring configuration to be stable. We see that the threshold 
electric field Ec decreases with increasing of the particle 
radius a. This tendency was experimentally observed in 
[38| for strong anchoring strength at the particle's surface 
that is in line with our speculations. 

But NLC with Ae < has the opposite tendency. If 
the particle has radius a > 0.16L and has Saturn ring 
director configuration the external field with value E > 

Ec = Et\J\ — { ^'^a^ Y makes hedgehog configuration to 
be stable. 

So it's possible to switch between Saturn ring and 
hedgehog configurations in the nematic cell with help of 
the electric (magnetic) field. Of course our speculations 
here are just a hypothesis and there should be experimen- 
tal or computer simulation test to check quantitatively 
the relation (|36l) . 
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FIG. 16. (Color online) (a) Log-log plots of the Fi as a function of the rescaled distance pjL. Dashed line 3 is the power-law 
asymptotics (b) Log-log plots of the D\ as a function of the rescaled distance pjL. Dashed line 3 is the power-law 

asymptotics |(^)'''. Planar cell, external electric field is applied along the y-axis, E|[t/ (see Fig[T]3). Blue lines 1 correspond to 
E — 0.99Et and Ae < 0. Green lines 2 correspond to the zero-field case E = 0. Red lines 4 correspond to E = 0.99Et and 
Ae > 0. This is a deconfinement effect in the planar cell. 




FIG. 17. (Color online) Map of the attraction and repulsion zones of interaction between two particles in the middle of the 
planar cell with rubbing along z and thickness L placed in the electric field parallel to no, E = 0.99i5t. Dashed lines are 
zone's borders in the unlimited NLC. Thick red lines a correspond to the positive Ae > 0, thick blue lines b correspond to the 
negative Ae < 0. Thin black lines are zone's borders in the zero-field case. Sign "-" means attraction, "-f " means repulsion, 
(a) Dipole-dipole interaction, pp' > 0. (b) Dipole-quadrupole interaction. The larger particle is in the center, pc > cp' . (c) 
Quadrupole-quadrupole interaction, cc' > 0. 



V. CONCLUSIONS 

In this paper we propose the Green function method 
which enables to describe quantitatively elastic colloidal 
interactions between axially symmetric particles confined 
in ncmatic cell under the action of the external electric 
(magnetic) fields. The current method is the generaliza- 
tion of the method proposed in p}3 . l33) for description of 
the colloids in confined NLC. General formulas for dipole- 
dipole, dipole-quadrupole and quadrupole-quadrupole in- 
teractions in the homeotropic and planar nematic cells 
with parallel and perpendicular field orientations are ob- 
tained. We consider cases of homeotropic as well as pla- 



nar nematic cells and have found some new results: 

1) Wc propose a criterion (|18p between radius of the 
particle and cell thickness when the Saturn ring configu- 
ration becomes stable ( for particle radius more than as) 
in the nematic cell in comparison with point hedgehog 
configuration and generalize this criterion on the case of 
the applied field ([55)) in the planar cell. Of course this cri- 
terion is not an exact result but rather an estimate which 
comes from analogy between exact results ([Ti)) - ([T5)) for 
unlimited NLC with external field and results for con- 
fined NLC with- and without external field. It should be 
tested with computer simulations or experiments. 

2) Deconfinement effect for dipole particles in 
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the homeotropic ncmatic cell with negative dielectric 
anisotropy Ae < and perpendicular to the cell elec- 
tric field when electric field is approaching it's Frederiks 
threshold value E => Ec- This means cancellation of the 
confinement effect for dipole particles (in [l3] the con- 
finement effect was found for quadrupole particles) near 
the Frederiks transition while it remains for quadrupole 
particles. The same effect takes place in the planar cell 
with Ae > 0; E_Lno and in the planar cell with Ae < 0, 
E||no. 

3) New effect of attraction and stabilization of the par- 
ticles in the homeotropic ncmatic cell along the electric 
field parallel to the cell planes when nematic dielectric 
anisotropy is negative Ae < 0. The minimun distance 
between two particles depends on the strength of the field 
and can be ordinary for Ae < (see Fig. [71). 

4) Attraction and repulsion zones for all elastic inter- 
actions are changed dramatically under the action of the 
external field. Especially when the electric field E is 
parallel to the planar cell planes and perpendicular to 
the rubbing direction no and NLC has negative dielec- 



tric anisotropy Ae < 0. The electric field widens lateral 
tails of the dipole-dipole interaction with increasing of 
the strength value E (sec Fig. [TO)) . But there is the 

critical value E^oi = 1.74£;t {Et = f y^jff ) at which a 

collapse of the lateral tales into dumbbell-shaped region 
occurs (compare Fig. [T^ and Fig. [TT]as well as Fig. [T3] 
and Fig. [H)) . Particles repel inside the dumbbell-shaped 
regions and attract outside of them (see Fig. [TT|) . So that 
the attraction in the perpendicular to no direction oc- 
curs at £■ > Ecoi ~ l.TAEt- Subsequent increasing of the 
field strength E compresses dumbbell-shaped regions and 
makes perpendicular attraction more significant. The 
same situation takes place for quadrupole-quadrupole in- 
teraction (see Figs. I13ll4ll5p . It is interesting that the 
same collapse and origin of the dumbbell-shaped regions 
occurs at the same field strength Ecoi = ^-T^Et as for the 
dipole-dipole interaction. 

We hope that listed above results may be tested either 
with help of computer simulations or experimentally and 
would be helpful in further investigations of fascinating 
world of colloids in liquid crystals. 
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